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We describe a lattice of asymmetrical qubit pairs in one or two dimensions, with couplings arranged 
so that the motion of single-qubit excited states mimics the behavior of charged lattice bosons 
hopping in a magnetic field. We show in particular that one can tune the parameters of the many- 
body circuit to reach a regime where the complex hopping phase between any two elements can be 
tuned by simply adjusting the relative phases of two applied oscillating voltage signals. We also 
propose a specific realization of our model using coupled three junction flux qubits, in which one can 
reach the strongly interacting bosonic quantum Hall limit where one will flnd anyonic excitations. 
The circuits could be used for topological quantum computation. 
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I. INTRODUCTION 

Fractional quantum Hall effects [H-O are among the 
most profound collections of phenomena to emerge in in- 
teracting quantum many-body systems. The elementary 
excitations in these systems do not act like bosons or 
fermions; rather, they are anyons, which in some cases 
can be used for a robust form of quantum computing 
[1, m . All physical examples of fractional quantum Hall 
effects are in two dimensional (2D) electron gases. Here 
we propose a method for linking standard qubit designs 
which will realize a bosonic fractional quantum Hall ef- 
fect. The rich theoretical literature on bosonic fractional 
quantum Hall effects suggests that there will be a large 
number of interesting states [6l-[l3| that could be explored 
in our system. These include 'Pfaffians' and their gener- 
alizations. Furthermore, one could anticipate that some 
important experiments (such as directly braiding quasi- 
particles) may be simpler in a qubit array than in a GaAs 
layer surrounded by AlGaAs. 

There are several competing approaches to engineer- 
ing bosonic fractional quantum Hall effects. One pro- 
posal uses Raman lasers to simulate the magnetic vec- 
tor potential in neutral cold atoms [13, [Hi- The tech- 
nical challenges are, however, quite daunting: new cool- 
ing methods need to be designed to offset heating from 
the Raman lasers, and the most natural probes are in- 
direct. Another scheme is to use lattices of tiny su- 
perconducting grains (charge qubits, [l6l - l2]| ) connected 
through Joscphson junctions. Suitably low temperatures 
can be reached in a dilution refrigerator, and the system 
is readily studied using transport measurements. Unfor- 
tunately, random charge noise, which scales linearly with 
the interaction strength, would prevent the quantum Hall 
regime from being reached without significant local tun- 
ing of the potentials on hundreds or thousands of lattice 
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sites. Other proposals include superconducting Jaynes- 
Cummings lattices j22| and "photon lattices" of coupled 
optical waveguides |23l - l25| , each of which have their own 
advantages and shortcomings. 

We here propose a new and promising approach. 
Namely, we consider a circuit of qubits, with a geometry 
which naturally maps onto a system of charged bosons 
hopping in a magnetic field. In order to produce complex 
hopping matrix elements we propose a lattice of coupled 
asymmetrical pairs of qubits, which we label as A or -B. 
We choose device parameters so that excitation energy 
lua of the A qubits is significantly smaller than the ex- 
citation energy of the B qubits, and place a B qubit 
on each link between neighboring A qubits. Further, we 
couple them to each other through alternating hopping 
(cr^cr^ +H.C., henceforth referred to as a "±" coupling) 
and potential {a^a'^g, a "zz" coupling) terms. We also 
apply an external oscillating electromagnetic field of fre- 
quency w to each qubit, with the relative phase of the 
signal applied to the B qubits shifted relative to that of 
the A qubits by a locally tunable ips ■ Since the B qubits 
arc higher energy than the A qubits, we can integrate 
them out, leading to complex tunneling matrix elements 
(the amplitude of a process where the states of neighbor- 
ing qubit pairs are exchanged) between A qubits with 
phases that can be tuned to any value by adjusting ifg. 

As we will describe below, a particularly attractive re- 
alization of this architecture would be to use three junc- 
tion "flux qubits" (FQs) (HH. The flux qubits are 
mesoscopic superconducting rings interrupted by three 
Joscphson junctions, placed in a magnetic field which is 
tuned so that nearly 1/2 of a magnetic fiux quantum pen- 
etrates the ring. The energies of the flux qubits can be 
tuned by adjusting this magnetic field, or by varying the 
areas of the Joscphson junctions, so that the B qubits are 
higher energy than the A qubits as outlined above. We 
then capacitively couple all the fiux qubits to an external, 
oscillating voltage Vq sinwt, and arrange the couplings so 
that the phases of the voltage applied to the B qubits 
are shifted relative to the A qubits. The subtle interplay 
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II. GENERAL FORMALISM 
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FIG. 1: Basic coupling structure for the A and B qubits. 
Each site in our many-body lattice would correspond to a 
single A qubit, which couples to its neighbors through one B 
qubit per link, joined through alternating hopping (±) and 
potential [zz) couplings as described in section [III Though 
drawn in one dimension in the figure, we ultimately intend to 
construct 2d lattices in this manner, and generalizations to 
even higher dimensions are also possible. 



of the oscillating applied voltages with the mix of charge 
(capacitive) and phase (Josephson) couplings introduces 
phase shifts which make these hopping matrix elements 
complex, mimicking the Peierls phases found for charged 
particles in magnetic fields. 

All of our flux qubits are operated in the regime where 
the Josephson energy i?j is large compared to the charg- 
ing energy Ec, so charge noise effects are exponentially 
suppressed. The system is therefore almost completely 
insensitive to stray low-frequency electric fields. The 
many-body excitation gap, a key feature of anyon states, 
can be measured through the single-qubit response to ap- 
plied oscillating voltages. The large nonlinearities of the 
flux qubit devices imply that the first excited states ex- 
perience an effectively infinite on-site repulsion. We note 
also that our scheme is not intended to function as a cir- 
cuit QED architecture (in contrast to the recent work 
of Koch et al [H, [s^] and others); the device parame- 
ters should be chosen so that the external voltages can 
be treated as purely classical sources, with no dynamical 
photons present in our system. 

The remainder of this paper is organized as follows. In 
section ini we write down the basic coupled qubit Hamil- 
tonian, and outline the conditions under which arbitrary 
external gauge fields can be simulated. In section Hill we 
describe three junction flux qubits, and how they can be 
coupled to obtain the arbitrary complex hopping phases 
derived in section [ill Having derived these phases, in 
section IIVI we show how the circuits of the two previ- 
ous sections can be used as building blocks for exotic 
boson fractional quantum Hall states. Finally, in sec- 
tion El we show how a simple arrangement of four flux 
qubits could experimentally demonstrate a nonzero effec- 
tive gauge field, and offer concluding remarks. 



A. Berry's Phase of a Rotating Spin 

Before outlining the physics of our qubit array, we 
would first like to discuss a simple example to more 
straightforwardly elucidate the origin of the complex hop- 
ping phases. Specifically, we will consider a pair of spins 
and examine the Berry's phase effects generated during 
a process where an excitation is transfered from one spin 
to its neighbor (whose eigenstates lie on a different axis 
from the first spin) by rotating both spins about z and 
then transferred back by rotating both spins about y. 
Specifically, let us consider two initially uncoupled spin- 
^ degrees of freedom, with the Hamiltonian, 

H ~ a\ + cos da% + sin 6 (cos Lpa^ + sin v^cr^) . (1) 

Let us assume that initially spin A is excited and spin B is 
in its ground state. We first act with the operator <j\(J% 
to transfer the excitation from Ato B, and assume energy 
is conserved in this process so that the final state after 
acting with (j\(Tg is |0a1b). We then act with (y\(Jg to 
transfer the excitation back; the resulting matrix element 
A4 for the entire process is 

M = (lAOsk>^|OAls)(OAlBkia||UOB) (2) 
= sin 6* (cos + i sin cos 61) . 

For B ^ njl and ^3 7^ 0, tt, is complex, and the re- 
sulting phase can be understood as a consequence of the 
Berry's phase acquired by a rotating spin, though we note 
of course that the Berry's phase discussed here is only an 
analogy, since we are considering the action of pairs of 
operators and not continuous, adiabatic changes to the 
system's wavefunction. When a spin m is rotated along 
a closed path, the resulting phase is equal to m times 
the area subtended by the path on the unit sphere. In 
this case, we have two spins which rotate, but both end 
in the same states in which they started, so we obtain 
a gauge-invariant phase equal to the sum of the phases 
picked up by both spins. The area subtended by A is 
just TT, but the area subtended by B depends on the pro- 
jection of Oy and onto its quantization axis, and thus 
depends on ips, yielding the result above. Note that if 
we'd acted with cf\a^Q or o\a\ twice instead of using a 
combination of the two, the outcome would necessarily 
be real, since M. would be the product of a matrix ele- 
ment and its Hermitian conjugate. In the Berry's phase 
picture, the phase is zero simply because the path of each 
spin's rotation would be a Id line, and thus each area is 
zero. Both inequivalent eigenstates and anisotropic op- 
erations are necessary for spin transfer matrix element to 
be complex. 

It is precisely this effect-the phase picked up by a spin 
which rotates as it propagates in space-which we will 
use to engineer artificial hopping phases in our lattice. 
Specifically, imagine the case in which we had two (iden- 
tical) A spins with a B spin in between them, and after 
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acting with cr^xcr^ to pass an excitation from Al to the 
B spin, we then act with o'^2'^b transfer the excitation 
to the second A qubit instead of sending it back to the 
first. Since the A spins are identical, the matrix element 
A4 should be the same as the one derived above, and 
therefore by letting B spins mediate a hopping coupling, 
we can introduce tunable phases in a lattice of A spins. 

Engineering this structure in a real spin (or qubit) lat- 
tice is by no means trivial. For real spins, one could in- 
troduce a spatially varying magnetic field to generate the 
inequivalent local eigenstates, but adding the anisotropic 
spin-spin interactions {cr^ag or cr\'J% instead of Sa • Sb) 
is very difficult. Conversely, for a more general lattice of 
qubits, generating passive anisotropic couplings is often 
straightforward, but generating inequivalent local eigen- 
states is not. Wc here demonstrate that coupling the 
qubits to a continuously oscillating monochromatic ex- 
ternal field can introduce the required rotations, provided 
that the phases of the signals applied to the B qubits are 
different from those applied to the A qubits. By adjust- 
ing these phases at a local level, we can independently 
tune the tunneling phase between any linked sites on the 
lattice, and can thus simulate any desired external gauge 
field, at least in principle. 

B. Qubit Coupling Hamiltonian 

We will consider a lattice of qubits, arranged such that 
there is an asymmetric pair of qubits A and B at each 
site. We shall assume throughout that the following con- 
ditions hold: 

(1) The nonlinearities of each physical system which we 
use as a qubit are large enough that we can consider them 
to be purely two-level systems, and ignore all eigenstates 
besides |0) and This requirement ultimately con- 
strains the magnitudes of the couplings between qubits, 
which must be small compared to the physical devices' 
absolute nonlinearities. 

(2) The qubits can be coupled to an external electro- 
magnetic field. We shall further require that the elec- 
tromagnetic field operator V (which could represent the 
coupling to magnetic fields as well) has no expectation 
value in either state, so (0| V |0) = (1| V |1) = 0. These 
fields will always be present in the qubit array Hamil- 
tonian, and wc will treat them in the standard rotating 
wave approximation. 

(3) We must be able to introduce two types of coupling 
between the qubits, so that the qubit-qubit Hamiltonian 
takes the form 

Hm = {a+<jg + <j^a+) + D''a%a%. (3) 

We must have independent control over both and 
for our method to succeed. Note that any physical 
coupling between the qubits will typically include terms 
which violate number conservation. However, when wc 
transform to the rotating frame when the external oscil- 
lating voltage is applied, the terms in Hi„t are unchanged 



but anomalous terms such as cr^cr;^ or crJ^cTg will be- 
come rapidly oscillating and can be dropped from the 
low-energy Hamiltonian. 

(4) We must be able to tune the relative phase ips of the 
external electromagnetic field applied to B qubits relative 
to the A qubits, as shown in fig.[TJ If 7^ 0, tt then time 
reversal symmetry is broken, since we cannot chose a zero 
point for the time t so that both Va (t) = Va {—t) and 
Vb {t) = Vb (— i). Breaking time reversal symmetry is a 
basic requirement for obtaining nontrivial effective gauge 
fields. 

These requirements could be fulfilled by a large number 
of physical systems, including spin qubits, trapped ions, 
and superconducting devices, which will be the focus of 
this work. Let us now consider the Hamiltonian of a 
given qubit pair. Has- Before turning on the oscillating 
fields, our qubit Hamiltonian is 

Hab = UJA(J\+^B(J% (4) 

Wc now turn on the oscillating fields. When acting on A 
or i?, we have: 

V ^ 2nA/B^%B^ (5) 

with 2^1a/b ~ {^a/b\ y |Oa/_b)i which we choose to be 
real (in the flux qubits described below, the full matrix 
element for the V operator is imaginary, but we have 
absorbed those factors of i into ct^). Wc now examine 

Vsinut = f^A/i3(e'"V;4/5+e-*"V+/5) (6) 
+nA/B [e a^i^ + e a^/^ j . 

Wc now transform to the rotating frame by applying the 
unitary transformation ji/;) — > exp— (cr;^ -|- cr^) i It/"). 
The time dependence of terms on the first line of © 
is cancelled out, leaving us with ^aib<^\ib P^^^ ^ 
of terms which are rapidly oscillating with frequency 2w. 
We now make the rotating wave approximation (RWA) 
to neglect these terms; if further accuracy is required we 
can treat them through a second order perturbation the- 
ory in and obtain a small correction to the 
'~'\lB terms. After transforming to the rotating frame 
and invoking the RWA, Hab is: 

Hab = {^A -oj)<Ja + {<^b - t^) (7% (7) 
+ ri/l(T^ -I- VIb {cosips<TB + s'mtpsCr^B) 
+ {D^ K'^B + or D^'a'X'^l,} . 

From now on wc will assume a; is tuned to resonance 
with the A qubits, so that w = lua and the single-site 
Hamiltonian for the A qubits is just flAO'A- 

To construct the full qubit lattice, we wire the qubits 
as in fig. [U so that the connection between any pair 
of neighboring A qubits consists of a zz coupling to 
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a B qubit followed by a ± coupling to the neighbor- 
ing A qubit. For simplicity, we will ignore cases where 
A qubits are coupled directly; such couplings will pro- 
duce either neighbor-neighbor potential interactions or 
real-valued hopping matrix elements, depending on their 
structure. We assume that the energy difference Eb — 

Ea ~ '2.\J {ujb — wa)^ + ^% — 2^ A = 5E is large com- 
pared to and -D^, so that wc can treat the A — B 
coupling perturbatively. We now eliminate the B qubits 
using second order perturbation theory; noting that all A 
qubits are identical, the resulting Hamiltonian, to order 
D'^ /5E, is given by: 

H = J2 {^^4^3 + H.c}j + 2^1a ^ ala„ (8) 

ij i 



COS! 



SE 

UJB - 



sin (cos fs{ij) + i cos 9 sin fs{ij) ) : 



^ [lJB - LOAf + 



Here aj/a,; creates/annihilates an excitation in the A 
qubit at site i, and VLa is equal to ^Ia plus O (J) shifts 
which depend on the coordination number of the lat- 
tice and magnitudes of the couplings. Since the qubits 
are spin-i, we have an effective hard-core constraint, so 

0. If we now identify 



aj II 



argj,., 



A • dr, 



(9) 



we see that the complex phases of J are identical to the 
Peierls phases of a charged particle moving on a lattice in 
an external gauge field A. Further, if we choose param- 
eters so that the B qubits are far off-resonance, will be 
small and 



oE 



(10) 



In this regime, we can freely adjust the phase of J without 
significantly altering its magnitude, and can thus simu- 
late any time-dependent external gauge field configura- 
tion we desire, simply by adjusting the B qubit phase 
shifts fsiij] each link. 

The ability to engineer artificial gauge fields of any de- 
sired configuration has tremendous potential to unlock 
new physics, and we will discuss the most natural ap- 
plication, simulating a uniform magnetic field to realize 
strongly interacting bosons in the quantum Hall regime, 
later in the work. Before doing so, however, we will first 
describe a possible implementation of this architecture in 
superconducting flux qubits. While fiux qubits arc cer- 
tainly not the only- or even necessarily the best- qubits 
to use for this purpose, our proposal will demonstrate 
that a fairly robust implementation of our architecture 
can be realized using device parameters from previous 
experiments. Thus, small lattices should be within reach 
of current technology. 



III. FLUX QUBITS 



The three-junction flux qubit consists of a supercon- 
ducting ring interrupted by three Josephson junctions as 
shown in fig. [5J with one junction whose area is rescalcd 
by a relative to the other two. A constant, tunable mag- 
netic flux bias of / 7^ 1/2 flux quanta is applied through 
the loop. Wc choose bottom third of the ring to be 
ground (which will be a physical ground in our case) with 
phase (/) = 0, then the two remaining degrees of freedom 
of the fiux qubit are the phases and (/)2 of the other 
two superconducting regions. The derivation of the fiux 
qubit Hamiltonian is descrbed in detail in Orlando et al 
[27| : in terms of the phases 0i and 02, the flux qubit 
Hamiltonian HpQ is 



Hfq — 



(1 + Q + 1]) {Ql + Ql) + 2aQiQ2 

{l + 7j){l + 2a + T])C 
—Ei [cos 01 -f cos (j)2 + a cos {2-kJ -f- 01 
2ti {aQ2 + {l + a + y-i)Qi) Vp sina;t 
(1 + r]){l + 2a + 77) 



(11) 



Here, Qj = —2eid/d<pj, Ej is the Josephson energy of 
the Josephson junctions and / is the total magnetic fiux 
through the loop in units of the magnetic flux quantum 
$0- The terms on the third line of ([TT|) represent the 
coupling of the flux qubit to the applied voltage Vq sinwt. 
For the moment, we will consider this Hamiltonian with 
1^0 = 0. 



We let 0± = (01 ± 02) /2. For / ^ 0, the symme- 
try between 0i and 02 is broken, and for / close to 
1/2, the ground and first excited states are distinguished 
by their behavior along the 0_ direction, as excitations 
along 0-|_ are significantly more expensive. The typical 
excitation energy for 0.4 < a < 0.6 and 0.5 </ < 0.55 
is ujpQ /2Tr = 12 - 30GHz for Ej/h - 200GHz and 
Ec ~ I2C = Ej/AO, and the nonlinearitics of the spec- 
trum arc all reasonably large. In this work we will only 
consider fiux qubits operated at the symmetry point of 
/ = 1/2, in which case the ground and first excited states 
are both even along 0+ and even or odd, respectively, 
along 0_. From this, we can readily translate operators 
in the phase basis to Pauli matrices acting in the qubit 
basis. We will define the following compact notation for 
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FIG. 2: (Color online) Basic circuit architecture. The regions enclosed in dashed boxes are three-junction flux qubits, which 
are connected to a physical ground. The blue {A, left and right) and red {B, center) qubits differ from each other by a rescaling 
of the area of the central Josephson junction, which is tuned so that the B qubits have higher energy excitations. A magnetic 
field penetrates the plane so that / flux quanta are enclosed by each ring. An oscillating voltage Ve (t) is applied near resonant 
transitions to both qubits, mixing their ground and first excited states. Excitations in the A fiux qubits can tunnel through 
the B qubits to each other; the oscillating voltage will make this transition matrix element complex. The qubit properties and 
the couplings between them are discussed in section Hill 



matrix elements; 



O.s 



{is\0\j,) e.g.A<'4',A= (OaIQiIIa). (12) 



In this notation, we have: 

Q, ^ 2e{-iyMllj7y 

sine/), ^ i-lYMZ^y^ 



(13) 



COSi 



COS 01 



^COS 01 ^ 



COS 01 



M 



00 

cos 01 



For consistency, all matrix elements Ai are calculated 
between the Vq = (non-rotating) cigenstates of the flux 
and transmon qubit Hamiltonians. 

Let us now turn to the coupling Hamiltonian between 
the qubits shown in fig. [2l We label the two A qubits 
by L and R. The coupling of the B qubit to the right 
qubit is a simple capacitive coupling, and so is given by 
a constant times ct^ct^, which becomes a ± coupling in 
the rotating frame: 



Hbr — 



9*, -B 



M 



(1 + 2a + 7?) (1 + 2/3 + 7?) 



.(14) 



It is important to note that both a^a^ and a^a^ become 
± couplings in the rotating frame (up to overall signs) , as 
the components of them which lead to net creation or de- 
struction of excitations arc rapidly oscillating and should 
be dropped. The coupling between the left qubit and 



the B qubit consists of two Josephson junctions; since 
these junctions define closed loops through ground, they 
pick up flux biases /' and /" from the external magnetic 
fleld. For simplicity, we choose the wiring geometry so 
that these biases are both zero mod 2tt. When we write 
the coupling between L and B as a set of Pauli matrices, 
the ± terms vanish due to the sign flips in p^ . but the 
zz term survives: 



H 



LB 



A/f™ "I 



(15) 



11 

cos 01 



M 



00 



3j (^L^l 



Alternately, one could obtain a pure zz coupling by sim- 
ply placing a single Josephson junction between a pair 
of regions, and choosing the wiring geometry so that the 
flux bias /' is nonzero, leading to an interaction term 
of the form —kEj cos {4>l2 — 4'Bi + 27r/') plus a capac- 
itive term with the same structure as p4)) . One could 
then tune /' so that the ± components of the xx and 
yy terms from these couplings interfere with each other, 
leaving only the zz part of the coupling. 

We are now in a position to plug in numbers and eval- 
uate J for this architecture. Consider flux qubits wired 
as in fig. [5J If we choose the realistic device parameters 
listed below, taking into account the single-qubit energy 
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|J|//i(GHz),argJ/7r 




FIG. 3: (Color online) Magnitude and phase of J for the 
device parameters given in eq. (|16|) at the resonance point 
u) = UJA- The blue and purple curves are |J| and aigj/ir, 
respectively, for Vb = 0.5mV; the yellow and green curves 
are the same quantities for Vb = ImV. J can be made sig- 
nificantly larger by increasing a and working away from the 
/ = 1/2 symmetry point, but the physical device nonlineari- 
ties are smaller in that regime and the system becomes more 
susceptible to fluctuations in the external magnetic field. 



shifts from the Z?^ coupling gives us: 



Ir. = 



Ei/h 



400nA, C = 3.25fF, a = 0.5, / = 0.5, (16) 

0.1, /3 = 0.45, K = 0.2, g = 0.2 

200GHz = 33Ec/h, uja^2tt x 15.5GHz, 

2^ X 18GHz, ^^^ = 2^x4.6^, 
Vo Vo mV 

2n X l.OGHz, 

27r X 1.4GHz. (17) 



A plot of J for Vb = 0.5, ImV is shown in figure fig. |3l 
calculated from For small values of Vb, \J\ is almost 
completely independent of ips , but for larger Vb the mag- 
nitude fluctuations become significant. \J\ can be further 
increased by up to an order of magnitude by choosing 
device parameters to work in the regime where / > 1/2, 
but the relative qubit nonlincarities arc smaller and the 
system becomes more susceptible to fluctuations in the 
external magnetic field. We emphasize that the parame- 
ters listed above certainly do not represent the best pos- 
sible choice for many-body physics, and indeed, it may 
ultimately turn out that other types of qubits may be su- 
perior for reaching the bosonic fractional quantum Hall 
regime described below. Nonetheless, they demonstrate 
that our system could be engineered with current tech- 
nology, and achieves hopping matrix elements which are 
around three orders of magnitude larger than the typical 
flux qubit decay and dephasing rates (around a MHz). 



IV. MANY-BODY STATES AND THE LOWEST 
LANDAU LEVEL 

By considering a 2d lattice of these qubits, requiring 
that I J| /{Eb — Ea) ^ 1, and ignoring inaccessible higher 



excited states, we arrive at the flnal hopping Hamiltonian 
dS]). Previous studies 0, i, [H!, have shown that 

the square lattice version of this Hamiltonian is analo- 
gous to the 2d lowest Landau level problem of strongly 
interacting bosons, and realizes abelian and non-abelian 
fractional ground states at the appropriate flxcd densi- 
ties. We expect that small arrays should be sufficient to 
observe quantum Hall physics, since the magnetic length 
Ib — 1/ V^TT^/ (where \I' is the gauge- invariant phase ac- 
cumulated when a particle circulates around a plaque- 
tte) can be less than a lattice spacing [HBl- Connec- 
tions between flux qubits beyond nearest neighbors can 
reproduce the exact lowest Landau level of the contin- 
uum [13, EHl and lead to more robust fractional quantum 
Hall states, but they may not be necessary to observe 
the Laughlin state at v = 1/2 Q- Here we adopt the 
standard definition of the filling fraction v as the ratio of 
particle to flux density. A wide range of other possible 
quantum spin-1/2 models with 2-body interactions, both 
with complex phases and without, could be studied in 
this device architecture; we find quantum Hall systems 
to be the most intriguing, due to the existence of abelian 
anyons at = 1/2 and the existence (with tuning) of non- 
abelian anyons ai v = \ and 3/2 [1,[4^, along with other 
exotic states at different filling fractions. The boson den- 
sity could be controlled by using a second external field 
at frequency lo' near 2U,a to populate the lattice; the lo' 
dependence of the system's response to this field could 
be used to measure the gaps of the many-body states. 

The incoherent particle gain and loss rates in our ar- 
ray from single qubit decay and dephasing effects should 
not be a significant obstacle to studying strongly cor- 
related many-body states. Using values from the pre- 
vious section and from the superconducting qubit lit- 
erature pll . |47| , a typical hopping parameter would be 
J/h = IGHz. It is important to note here that the 
particles in our case are rotating frame antisymmetric 
superpositions of the ground and first excited states of 
the physical qubits, so their decay would be from one 
type of superposition to the other. This rate would be 
roughly given by the dephasing rate of the qubits, which 
for flux qubits is of order IMHz. With a Landau band 
spacing of ojlll — 3J in a square lattice at = 1/4 
quanta per plaquette, the relative correction to the Lan- 
dau bandwidth from this process decay would be thus be 
insigniflcant. We expect that this loss rate by itself will 
not prevent quantum Hall states from forming in our ar- 
ray [2^ . Likewise, a small number of "dead" sites (where 
a qubit is defective and cannot be excited) should also 
be relatively harmless- the many-body wavcfunction can 
eliminate these defects simply by nucleating a quasihole 
at each site. So long as the density of flux quanta is 
large compared to the defect density, these defects will 
simply make small shifts in the gap energy and particle 
density of the gapped states, but will not have any other 
qualitative effects on the system. 

More worrisome is the issue of time-independent ran- 
dom variations in the qubit properties at every site, which 
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could disrupt the formation of topological states if these 
variations became large enough. To quantify this issue, 
we numerically simulated the broadening of the lowest 
Landau level in our model as a function of three static 
(quenched) noise sources: random fluctuations in the on- 
site potential (shifts in the rotating frame excitation en- 
ergy of a given A qubit), random fluctuations in the mag- 
nitude of J, and random fluctuations in the phase of J 
between neighboring sites. In a real system, these noise 
sources would be correlated, but as the details of those 
correlations would depend in part on the physical imple- 
mentation of the qubits, we have assumed that each type 
of quenched disorder is applied randomly to every site 
with no dependence on the other types or on the disor- 
der at nearby sites. To determine the broadening from 
each noise source, we numerically diagonalizcd the single- 
particle hopping matrix on 8x8 and 12x12 lattices with 
periodic boundary conditions, given by the Hamiltonian: 

Hlll = (e'(^"+"'*-) + H.c}j (18) 



E 



JNNSUiTli 



Here, the hopping matrix elements are restricted to near- 
est and next nearest neighbors with relative magnitudes 
chosen as in j^], 5Ui and 54>ij are dimensionless param- 
eters which are Gaussian distributed about 0, Jnn is the 
average nearest neighbor hopping energy, and Fij is a 
dimensionless parameter which is Gaussian distributed 
about 1. We diagonalized ([TS]) for 25 random distribu- 
tions of noise for each data point, and from the spec- 
trum we extracted the lowest Landau level broadening 
A, which is the ratio of the energy splitting between the 
lowest and highest LLL states to the splitting between 
the highest LLL state and the bottom of the first ex- 
cited band. We then fit A [cfu/j/^) as a function of the 
standard deviation of each noise source with the other 
two sources set to zero; this relationship was linear in 
each case for small fluctuations. The results of our sim- 
ulations are shown in table H] note that Aq is nonzero 
even without defects, as a consequence of truncating the 
Hamiltonian in [43 | to nearest and next nearest neighbor 
hopping. 

It is important to note that this calculation only cap- 
tures distortions to the single particle spectrum and that 
the many-body response to noise of this type is a sub- 
tle problem beyond the scope of this work. However, we 
qualitatively expect that the topological states should 
be disrupted when the normalized Landau level splitting 
A approaches the dimensionless quasiparticle excitation 
gap Aqp/ Jnn . In numerical studies of this system in 
the clean limit with hard-core 2-body repulsion (largely 
unpublished), Aqp/J^N typically ranged between 0.2 
and 1 for correlated states at different flux and parti- 
cle densities, and tended to be larger at smaller filling 
fractions. This suggests that many-body quantum Hall 
states should exist in our system when noise is sufficiently 
well-controlled. 



Flux Density 


Ao 


Cu 


Cj 




1/4 


0.015 


0.41 


1.42 


1.75 


1/3 


0.018 


0.72 


1.21 


2.36 


3/8 


0.08 


0.35 


0.99 


1.94 



TABLE I: Robustness of the lowest Landau level to external 
noise sources. For the random noise simulations described in 
the text, we fit the normalized splitting A of the lowest Lan- 
dau level to the function A = Aq -I- Cucru + Cjaj + C*cr*, 
where the cr's are the standard deviation of each noise source 
(local potential, hopping magnitude, and hopping phase) 
which is applied randomly to every site (au) and link be- 
tween sites (crj and cr*). As seen in the Hamiltonian p8[) . 
the potential fluctuations are in units of Jjvjv and the phase 
fluctuations are in units of tt. Above the flux density ^' = 1/3, 
truncation to nearest and next nearest neighbor hopping in- 
troduces significant broadening even in the clean system, so 
flux densities of 1/3 or less should be the focus of experiments 
on our design. 



V. A SIMPLE EXPERIMENT TO 
DEMONSTRATE THE GAUGE FIELD 

While the ultimate purpose of this proposal is to study 
exotic many-body states in an array of hundreds or thou- 
sands of flux qubits, the existence of a nontrivial gauge 
field can be demonstrated by studying an arrangement of 
four fiux qubits, connected in a loop. Consider a square 
loop of four flux qubits labeled (1-4), where qubit 1 sits 
at the top left corner and qubit 4 at the bottom right, 
as shown in Fig. U) For this choice, any hop through a 
coupling will accumulate a phase ip, giving a total of 

= 2ijj for a complete circuit of the loop. Conversely, 
if the phases of the voltages applied to the B qubits are 
shifted by tt from one FQ-FQ pair to the next, the mag- 
nitude of the hopping matrix element will be unchanged 
but there will be no complex phase accumulation. In this 
case, the B qubits have identical rotating frame energies 
to the A qubits, and differ from them through the rela- 
tive phases ipsi of the applied voltages. We will assume 
for simplicity that the magnitudes of the hopping matrix 
elements from the D'^ and couplings are both equal 
to J. 

To demonstrate that the alternating voltages generate 
a nonzero effective fiux through the four-qubit loop, we 
first initialize the array by letting all four qubits relax to 
their ground states. At time < = 0, we apply a microwave 
pulse to qubit 1 to excite it into the fluxon state 1 1) , and 
then at time t we measure the state of qubit 4. The 
probability of qubit 4 being occupied by the fiuxon is 
given by 



Piit) = (OiOaOaUle 



iHt/h 



I1O2O3O4) 



(19) 



1 



2tJ * 

, cos I COS — 

4 V \ h 4 



f2tJ . 

cos sm — 

\ h 4 



This interference pattern is particularly striking when ^ 
is nearly equal to tt. If we let ^' = tt -|- e, the probability 
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distribution becomes 



Pi (t) = sin 



V2Jt \' 



Jte 



sm ■ 



-V 



fi J \ 2V2hJ 



(20) 



In the limit of e — > 0, the probability of qubit 4 being 
occupied becomes zero at all times, due to the perfect 
interference of the two paths. This is a dramatic effect, 
and while field fluctuations and fabrication defects would 
prevent perfect interference in a real device, the strong 
slowing of the occupation periodicity of qubit 4 as 4* ap- 
proaches TT would be readily observable. Such interfer- 
ence is only possible if there is a gauge-invariant phase 
difference between the two paths, and would therefore 
demonstrate that nontrivial effective gauge fields are re- 
alized in our architecture. 



VI. CONCLUSION 

We have demonstrated a method for realizing a quan- 
tum Hall state of bosons using asymmetric qubit pairs, 
driven by applied oscillating electric fields. We also 
demonstrated that our model could be implemented in 
a lattice of flux qubits. With appropriate protocols for 
stabilizing the average particle density and measuring the 
conductivity, we expect that conductivity quantization 
could be observed on small arrays, though we note that 
the details of how to measure the conductivity are beyond 
the scope of this article. The statistics of anyonic collec- 
tive modes could be determined through similar methods 

[Mill. 

Further, the dynamical tunability of our model could 
be exploited to realize exotic combinations of states that 
would be difficult or impossible to study in cold atom 
or solid state systems. One could locally adjust the ap- 
plied voltage Vq sin {ojt + if s) and flux bias / to change 
the gauge field density and effective chemical potential in 
a given region, creating islands of arbitrary shape which 
could be at a different filling fraction than the surround- 
ing lattice and thus have different anyonic modes. Alter- 
nately, by reversing the signs of all the phase shifts ips 
in a region, one can crate a sharp boundary between re- 
gions with effective gauge fields of equal magnitude but 
opposite sign. In both cases, we expect physics along the 
boundaries to be rich. 

Finally, by locally tuning Vq, (fs and / to manipulate 
vortices in the qubit lattice, arrays of ordinary qubits 
could be used to construct a topological non-abelian 
anyon qubit 0, S Is^ . trading information density for 
topological protection against decoherence. In that sense 
our proposal is similar to the surface code and cluster 



state [53|, |54| ideas developed in recent years, and pro- 
vides a new potential mechanism for reducing decoher- 
ence in superconducting quantum information devices. 



Vb sin ojt 







^ 



Vb sin (ujt + ifs) 



gc 



HI" 











FIG. 4: (Color online) Configuration to demonstrate the arti- 
ficial gauge field, as outlined in section|Vl As described in sec- 
tion[Vl appropriately tuning the phase ofi'sets (fsi will produce 
a gauge- invariant phase difi'erence in the two paths that the 
mobile fiuxon excitation could take from qubit 1 to qubit 4. 
The resulting interference of these two paths can be detected 
in the time-dependent probability P4 {t) of qubit 4 being in 
its excited state. 
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